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Анотація. Діоксид сірки - один з газів, що найчастіше зустрічається, який забруднює повітря та несе 
велику кількість небезпек для здоров'я людини і довкілля. З метою зниження шкідливої дії, необхідно 
контролювати викиди димових газів на електростанціях, оскільки значна частина діоксиду сірки в атмосфері 
виникає з димових газів, що утворюються в процесі виробництва електроенергії. У роботі розглядаються 
проблеми оптимізації процесу десульфуризації димових газів від діоксиду сірки за допомогою інтелектуального 
аналізу даних. Розроблено модифікований генетичний метод оптимізації процесу десульфуризації димових газів 
на основі нейронної мережі, що дозволяє відображати вплив характеристик часового ряду параметрів процесу на 
фактичну ефективність десульфуризації 1 підвищувати її точність прогнозу. Принципова відмінність 
запропонованого генетичного методу від існуючих аналогів полягає у використанні адаптивної мутації, яка 
використовує в процесі своєї роботи рівень розвитку популяції. У ній найменш значущі гени з більшою 


59 


155 2710 - 1673. Агбіїсіа! Паліеїйдепсе. 2021. Мо 1 


ймовірністю будуть мутувати в хромосомах, ніж з високою придатністю, це підвищує точність та збільшує їх 
роль в пошуку. Було проведено порівняльний аналіз розробленого генетичного методу та відомих, в результаті 
чого було встановлено, що він дає найменшу помилку прогнозування вмісту 50» на виході та дозволяє зменшити 
час виконання прогнозування ефективності десульфуризації димових газів. Отримані результати дозволяють 
запропонувати дійовий метод для підвищення ефективності процесу десульфуризації димових газів та зниження 


викидів 50» в атмосферу. 


Ключові слова: знесірчення димових газів, діоксид сірки, штучна нейронна мережа, генетичний 


алгоритм. 


Вступ 

За | даними | Майше  Сеобзсіепсе, 
супутники МАЗА виявили 500 нових 
джерел забруднення повітря, близько 40 з 
них небезпечні діоксидом сірки |1|. Ця 
речовина вважається одним з найбільш 
ризикованих для атмосфери Землі. 
Основними джерелами викидів двоокису 
сірки | є теплові електростанції, що 
працюють на твердому й рідкому паливі, а 
також металургійні підприємства. Отже, 
контроль вмісту 502, в димових газах від 
спалювання вугілля є ефективним методом 
щодо скорочення викидів в атмосферу |21. 

Викиди 50» можна зменшити за 
рахунок встановлення обладнання для 
сіркоочищення |на новозбудованих та 
існуючих вугільних блоках, а також 
відповідних вимог до ефективності 
знесірчення. Щоб підвищити ефективність 
знесірчення і більш ефективно знизити 
викиди 502, | необхідно | додатково 
оптимізувати систему управління 
знесірчення, відповідно до потреб галузі 
Г21. 

Точне встановлення взаємозв'язку 
між параметрами процесу і фактичною 
ефективністю сіркоочищення є основою 
для оптимізації з системи управління 
сіркоочищенням. В даний час велика 
кількість локальних приладів для аналізу і 
моніторингу | димових газів | може 
безпосередньо контролювати 
концентрацію 502 в димових газах на вході 
й виході обладнання для десульфуризації 1 
розраховувати її ефективність. Однак цей 
метод являє собою всього лиш простий 
зворотний зв'язок за результатами процесу 
реакції знесірчення 1 не відображає 
результатів моніторингу. У той же час на 
обладнання для моніторингу сильно 
впливають зовнішні фактори, іноді 
виникають збої і тому подібне, що 
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призводить до неточних результатів 
вимірювань. Отже, на ефективність 
знесірчення впливають різні фактори, такі 
як значення рН суспензії, температура 
димових газів на вході, щільність суспензії 
абсорбційної вежі і масова концентрація 
50» на вході димових газів |21. 

За останні роки були запропоновані 
численні підходи, включаючи 
експериментальні дослідження, 
математичні моделі та моделі машинного 
навчання, для прогнозування ефективності 
десульфуризації. Серед них - математичні 
моделі 1 моделі машинного навчання, а 
також їх комбінації викликали великий 
дослідницький інтерес. Однак 
прогнозування ефективності 
десульфуризації все ж складно повністю 


змоделювати математично. Деякі 
дослідження спрощують систему за 
допомогою | припущень, які вносять 


помилки в прогнози. Крім того, чисельні 
рішення, які використовуються в цих 
математичних моделях, вимагають великих 
обчислювальних ресурсів. | У цьому 
контексті моделі машинного навчання 
виявляються корисними інструментами для 
прогнозування |21. 

Один з перспективних напрямків 
вирішення даної проблеми трунтується на 
застосуванні штучних нейронних мереж 1 
генетичних алгоритмів, як найбільш 
прогресивних по відношенню проблем 
прогнозування ефективності 
десульфуризації димових газів. Генетичні 
алгоритми відносяться до класу методів 
пошуку, які ітераційно поліпшують якість 
варіантів вирішення процедурами 
рекомбінації 1 відбором для виживання. 
Завдяки універсальності обчислювальної 
схеми, можливостям паралельної реалізації 
та стійкості до шуму, ГА знаходять успішне 
практичне застосування при вирішенні 
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багатьох складних нелінійних 
багатовимірних задач оптимізації | 3). 

У | цій | роботі | запропоновано 
модифікований генетичний метод 
оптимізації | процесу | десульфуризації 
димових газів на основі нейронної мережі, 
що дозволяє відображати вплив 
характеристик часового ряду параметрів 
процесу | на | фактичну ефективність 
десульфуризації і підвищувати її точність 
прогнозу. 


Постановка завдання 

У наш час найбільший інтерес являє 
використання технологій десульфуризації 
пилогазоповітряних | потоків, | здатних 
забезпечити ефективне зв'язування 
діоксиду сірки. Одним з найбільш 
поширених методів є метод мокрого 
сіркоочищення |4|. В основі методу 
мокрого абсорбційного очищення димових 
газів лежить нейтралізація | сірчистої 
кислоти, що виходить в результаті 
розчинення діоксиду сірки, що міститься в 
димових газах, гідратом оксиду кальцію 
Са( ОН)» (вапном) або карбонатом кальцію 
СаСО.  (вапняком). Вапно | - більш 
дорожчий і дефіцитний матеріал, однак 
застосування вапна в якості сорбенту 
дозволяє майже в два рази скоротити 
витрати | реагенту, знизити витрати 
електроенергії на приготування суспензії 1 
зрошення абсорбера. Ефективність 
використання Са(ОН)» для очищення газу 
залежить від температури газу, Його 
вологості й часу перебування реагенту в 
системі. Для регулювання температури 
застосовується інжекція дрібнорозпиленої 
води | під | високим | тиском, яка 
випаровується | 1 | прискорює реакції 
абсорбції |51. 

Проаналізувавши принцип роботи 
методу мокрої абсорбційної 0 очистки 
димових газів, було встановлено, що 
узагальнену модель залежності ступеня 
сіркоочищення від використовуваних 
реагентів з деяким припущенням, можна 
призвести до виду (1) | 6: 
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К Фе5 - і ( Хр - топі а45? Хр -абег о а45? Х Лі маре? 


, (1) 
фена 


Х опі зіаре? Хши 9" гесіг У "умер? Хіте 7 Хо оп геарепі ) 


ДО | Хегр йопг аз ї - | ТЄмпература перед 
абсорбером, Хстр. айег аз - Температура після 
абсорбера, Ха зчіазг - Об'єм газу на 1-у 
СХОДИНКУ, Хкесопа зіаєе - об'єм газу на 2-у 
СХОДИНКУ, Хавесопаі 5іаре - запиленість після 2-Ї 
СХОДИНКИ, | Хиесіг об'єм газу | на 
рециркуляцію, Хуаіег витрати води у 
абсорбері, Хіпе- ВИТрати вапна, Хуєсопа геавеп - 
витрати вторинного реагенту. 

Отже, для вирішення даного завдання 
необхідно провести аналіз результатів 
роботи  десульфуризації димових газів, 
встановити залежність | між ступенем 
сіркоочищення від обсягу 
використовуваних реагентів та створити 
модель цієї залежності, яка надала б 
можливість прогнозувати ефективність 
знесірчення димових газів, внаслідок чого б 
знизилися викиди 502 та експлуатаційні 
витрати | за | рахунок оптимального 
використання реагентів у майбутньому. 


Аналіз літературних 

постановка проблеми 

У статті |7| пропонується модель 
прогнозування ефективності  десульфу- 
ризації димових газів, заснована на 
нейронній мережі І.5ТМ. Автори узяли в 
якості розгляду 1000МУУ установку в Китаї, 
розглядаючи | основні фактори, що 
впливають на ефективність  десульфу- 
ризації вологих димових газів вапняком 1 
гіпсом. Для побудови моделі 
прогнозування ефективності  десульфу- 
ризації використовується нейронна мережа 
з короткочасною пам'яттю (1.5ТМ). 

Загальна структура моделі 
прогнозування І5ТМ складається з п'ять 
функціональних модулів: вхідний рівень, 
прихований рівень, вихідний рівень, 
навчання мережі та прогнозування мережі. 
Вхідний рівень відповідає за первинну 
обробку | вихідного | часового | ряду 
несправностей, щоб задовольнити вхідним 
вимогам мережі. 


даних та 
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Прихований рівень використовує 
комірки І.5ТМ для побудови одношарової 
циклічної нейронної мережі, вихідний 
рівень надає результати прогнозування. 
При навчанні нейронної мережі ІЗТМ 
використовується з алгоритм зворотного 
поширення за нчасом  (Васкргорабайоп 
іФгоцеб їйте), який аналогічний класичному 
алгоритму зворотного поширення. 

Перевагою запропонованої моделі в 
порівнянні з іншими моделями є те, що 
вона може відображати вплив 
характеристик часових рядів параметрів 
процесу | на | фактичну ефективність 
десульфуризації і підвищувати точність 
прогнозування. 

Недоліком запропонованої моделі є 
те, що під час навчання складно трактувати 
результати, це накладає обмеження на 
можливість поліпшення моделі, так само 
неможливо передбачити, коли зміниться 


динаміка |  десульфуризації | 1 | модель 
закінчить працювати. 

У роботі 131 представлено 
прогнозування ефективності мокрої 


десульфуризації вапняку-гіпсу на основі 
методу опорних векторів (5УМ). Автори 
роботи приймають співвідношення рідина- 
газ, швидкість димового газу, об'єм 
окислювального повітря, температуру газу, 
пил, що міститься в димовому газі, 
концентрацію діоксиду сірки на вході, рН 
суспензії абсорбційної вежі, концентрацію 
суспензії | абсорбційної | вежі,  спів- 
відношення кальцію і сірки як незалежні 
змінні, а ефективність десульфуризації - як 
залежну змінну. Ідея запропонованого 
методу полягає в побудові гіперплощини, 
яка виступає в якості поверхні рішень, 
максимально розділяє позитивні і негативні 
приклади з навчальної множини. Машина 
опорних векторів є  апроксимуючою 
реалізацією методу мінімізації структур- 
ного ризику, який заснований на тому, що 
рівень помилок навченої машини на 
тестовій множині представлений у вигляді 
суми помилки навчання та доданка, який 
залежить від вимірювання показника 
Вапніка-Червоненкіса. 
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Перевага | розглянутого методу 
полягає в тому, що для вирішення задачі 
класифікації, на відміну від більшості 
інших методів, для нього достатньо 
невеликого набору даних. 

Недолік методу полягає в тому, що 
для класифікації використовується не вся 
множина зразків, а лише їх невелика 
частина, яка знаходиться на кордонах. 

У роботі 191 представлено 
дослідження з інтелектуального аналізу 
даних 1 оптимізації роботи систем мокрого 
сіркоочищення димових газів. Авторами 
був запропонований комплексний критерій 
оцінки ї з використанням мінімальних 
витрат в якості цільової функції для 
отримання оптимальних умов роботи 
системи мокрого сіркоочищення димових 
газів. Для підвищення точності 
інтелектуального аналізу цільових даних 
був запропонований поліпшений метод 
нечіткої кластеризації. Цей алгоритм 
використовував результати роботи методу 
К-середніх як початкових умов та методу 
нечітких  С-середніх як аналітичного 
методу. Процедура вирішення  запро- 
понованого методу виражена у два етапи: 
використання методу  К-середніх 0 для 
отримання початкової кількості кластерів 
та їх центрів, а також використання 
нечітких  С-середніх для розрахунку 
остаточних результатів кластеризації. 

Перевагами запропонованого методу 
є відносна простота реалізації, відсутність 
необхідності експертної розмітки 
навчальної вибірки і відносно невисокі 
вимоги до обчислювальних потужностей. 

Недоліком запропонованого методу є 
його низька стійкість до викидів (у даних, 
які можуть спотворювати середнє). Як 
варіант, вирішенням цієї проблеми може 
бути використання модифікації алгоритму 


К-медіани. Також недоліком є те, що 
результат його роботи залежить від 
початкового вибору  центроїдів і у 
загальному випадку метод не дає 


оптимальних результатів, а дозволяє знайти 
тільки субоптимальні рішення. 
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У роботі |10| представлено вирі- 
шення завдання скорочення викидів для 
системи примусового окислення димових 
газів  вапняком | на основі штучного 
інтелекту. Авторами була створена штучна 
нейронна мережа на основі багатошарового 
перцептрона (МІР) для дослідження 
взаємозв'язків між вхідними керуючими 
змінними | та | вихідними змінними. 
Запропонована модель продемонструвала 
здатність вивчати складну кінетику реакції 
між реагентами. Вхідними параметрами 
моделі є: рН суспензії абсорбера, маса 502 
на вході, температура, маса МО» на вході, 
відсотковий вміст кисню, обсяг повітря для 
окислення, щільність суспензії абсорбера, 
об'ємний вміст вологи, маса часток пилу як 
доля об'ємної витрати димових газів. 
Вихідними параметрами моделі є: маса 
50», Не, МО» і пилу на виході. Модель 
складається з одного прихованого шару, 
який містить 27 нейронів. Навчання 
відбувається | за | допомогою методу 
зворотного поширення помилки. Воно 
виконується до тих пір, поки не буде 
виконано один із критеріїв зупинки, тобто 
буде досягнута зміна помилки збіжності на 
0,0000001, або під час навчання буде 
досягнута максимальна кількість епох, за 
межами яких відбувається незначне 
зменшення помилки. 

Перевагою запропонованого методу є 
те, що він має високий ступінь зв'язності, 
який | реалізовується | за | допомогою 
синаптичних з'єднань. Також при навчанні 
представленої нейронної мережі, за 
рахунок свої внутрішньої будови, вона 
встановлює закономірності у зв'язку 
вхідних і вихідних даних, тим самим як би 
"узагальнює" отриманий на навчальній 
вибірці досвід. 

Недоліком запропонованого методу є 
те, що він використовує зворотнє 
поширення помилки, в результаті чого 
процес навчання триває досить довго. 
Надмірне збільшення ваг нейронної мережі 
може призвести до її «паралічу», коли 
синаптичні ваги перестануть змінюватися, 
не досягнувши кінця навчання. 
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Величину кроку необхідно постійно 
підлаштовувати, мережу необхідно 
постійно  «донавчати», так | як вона 
поступово «забуває» попередні навчальні 
набори, якими вона була навчена. 

У роботі (11) представлено вирішення 


задачі | оптимізація | системи мокрої 
десульфуризації димових газів. | Для 
прогнозування | викидів | 502, була 


розроблена модифікована модель мокрого 
сіркоочищення димових газів, що поєднує 
математичну модель та штучну нейронну 
мережу. 

Було запропоновано модифікований 
Р5О  (Рагіісіє  5мапаі  Оріїтігайоп) з 
алгоритмом | штрафної функції для 
підвищення економічності процесу 
десульфуризації, спрямованого на 
мінімізацію | експлуатаційних витрат. 
Метод | рою часток заснований на 
обчисленні в кожній ітерації координат 
центроїда | 1 | швидкостей. | По суті, 
вирішується завдання мінімізації фітнес- 
функції: чим менше значення фітнес- 
функції на кожній ітерації, тим ближче 
результат до оптимального центроїду рою 
частинок. Фітнес-функція визначається як 
середнє значення евклідових відстаней між 
вектором частки та вектором, що визначає 
координати центроїда. 

У результаті чого були оптимізовані 
робочі | параметри, включаючи рН 
вапняного розчину, температуру, щільність 
вапняного розчину і кількість працюючих 
циркуляційних насосів. Результати 
показують, що Р50О з функцією штрафу 
забезпечує задовільну продуктивність 
збіжності для зниження експлуатаційних 
витрат при обмеженнях. 

Недолік даного методу полягають в 
тому, що він легко потрапляє в локальний 
оптимум в багатовимірному просторі та 
вимагає великі часові витрати на пошук 


рішення. 

У роботі П12| представлено 
прогнозування та оптимізацію роботи 
системи десульфуризації з використанням 
нейронної | мережі та генетичного 
алгоритму. 
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Автори роботи, приймаючи 
коефіцієнт десульфуризації та економічну 
вартість в якості двох цілей, створили 
модель прогнозування з десятьма входами і 
двома виходами. Був створений 
модифікований генетичний алгоритм для 
побудови моделі та оптимізації вартості 
процесу  десульфуризації. Генетичний 
алгоритм складається з двох етапів. 

На першому етапі ГА знаходить 
субоптимальну архітектуру ШНМ. На 
другому етапі визначаються субоптимальні 
значення вагових коефіцієнтів та зсувів 
ШНМ. Перший етап можна вважати 
підготовчою стадією навчання, завдяки 
якому спрощується подальше завдання 
визначення ваг і зміщень. Авторами був 
використаний наступний підхід: хромосома 
будувалася з числа нейронів у прихованому 
шарі, імпульсу, швидкості навчання, типу 
функції активації й алгоритму навчання 
ШНМ та декількох інших параметрів. У 
зв'язку з різноманітністю даних параметрів, 
використовувати оператор інверсії в них 
неможливо в загальному випадку. У число 
додаткових генів був включений вектор, що 
характеризує входження тієї чи іншої 
змінної в навчальний масив. Нулі й одиниці 
розподілялися по локусам випадковим 
чином, з тим лише зауваженням, що число 
нулів не повинно перевищувати половину 
вихідного числа змінних. Введення такого 
типу генів у хромосому проводилося з 
метою виявлення субоптимального набору 
змінних. Після визначення архітектури 
ШНМ здійснюється 0 запуск навчання 
мережі. Після чого формується ряд векторів 
ваг і зміщень, приведених до розміру пхіІ, і 
вектор  середньоквадратичної помилки 
навчання. З ряду векторів (хромосом) 
створюється початкова популяція, яка 
згодом буде еволюціонувати за допомогою 
виконання генетичного алгоритму. 

Перевагою запропонованого методу є 
підвищення точності класифікації при 
використанні ГА-модифікованого підходу. 
ГА дозволяє досвідченим шляхом підібрати 
найкраще поєднання основних ознак 
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ШНМ, яке може варіюватися не тільки в 
різних / завданнях, але й при зміні 
статистичних параметрів даних. Таким 
чином, даний підхід дозволяє оптимізувати 
процес десульфуризації. 

З недоліків даного методу можна 
виділити витрачання набагато більшого 
обсяг обчислювальних ресурсів. Тому 
представлений ГА можна використовувати 
тільки | при наявності | доступу до 
розподілених обчислювальних ресурсів. 

Аналіз | робіт  |7-12| дозволяє 
стверджувати, що проведення досліджень 
щодо оптимізації та прогнозування 
ефективності очищення димових газів від 
оксидів азоту є доволі актуальною задачею. 
Також у ході роботи було встановлено, що 
найбільш ефективним методом для 
вирішення поставленого завдання є 
використання нейронних мереж. Оскільки 
вони мають такі переваги: можливість 
адаптації до змін, яка дозволяє Їм 
працювати навіть у критичній ситуації, 
самонавчання, швидкість роботи 
(складання прогнозу й планування може 
йти набагато швидше, ніж при 
використанні точних алгоритмів), а також 
висока стійкість до шумових даних та 
можливість застосування необмеженої 
кількості незалежних змінних. При цьому, 
незважаючи на їх очевидні переваги, 
основним недоліком є складність вибору 
ефективної | методики | навчання та 
створення навчального масиву даних, а 
також складність правильного розподілу 
ваги матриці на початковому етапі, в 
умовах невизначеності. Для вирішення цієї 
проблеми добре підійде генетичний метод. 
Перевага даного методу полягає в більшій 
точності та зниженні кількості часу, що 
витрачається на навчання. Він ефективний 
для навчання нейронної мережі типу 
«багатошаровий перцептрон», «ШНМ (із 
загальною регресією», «мереж Кохонена» 
та ін. Застосування еволюційних гене- 
тичних алгоритмів на етапі навчання 
допомагає скоротити кількість навчальних 
циклів за | рахунок використання 
параметрів, що мають найбільшу вагу. 
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Отже, з метою вирішення проблем, 
що | виникають | при використанні 
досліджуваних методів, було прийнято 
рішення розробити модифікований 
генетичний метод для налаштування 
вагових коефіцієнтів штучної нейронної 
мережі для підвищення ефективності 
процесу десульфуризації димових газів та 
зниження викидів 502 в атмосферу. 


Мета й задачі дослідження 


Об'єкт дослідження  - системи 
десульфуризації на вугільних 
електростанціях. 

Предмет дослідження  - методи 
прогнозування | ефективності | десуль- 
фуризації димових газів. 

Методи дослідження - традиційні 


моделі (дерева рішень, метод найближчого 
сусіда, метод мурашиної колонії), нейронні 
мережі, комбіновані методи (нейронні 
мережі й генетичні алгоритми, нейронні 
мережі й мультиагентні системи). 

Мета дослідження | - розробка 
математичної моделі вирішення задавання 
оптимізації | процесу | десульфуризації 
димових газів від діоксиду сірки з 
використанням нейронної мережі. 

Для досягнення мети дослідження 
необхідно було розв'язати такі завдання: 

- розробити модифікований  гене- 
тичний метод навчання нейронної 
мережі з використанням адаптивної 
мутації для прогнозування 
ефективності знесірчення димових 
газів; 

- представити експериментальне 
дослідження запропонованого 
генетичного методу. 


Розробка генетичного методу для 


оптимізації процесу  десульфу- 
ризації димових газів від діоксиду 
сірки 


Найбільш відомим з еволюційних 
методів в теперішній час є генетичний 
алгоритм для знаходження глобального 
екстремуму багатоекстремальної функції. 
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Він полягає в паралельній обробці 
множини альтернативних рішень, при 
цьому пошук концентрується на найбільш 
перспективних з них. Це свідчить про 
можливість використання генетичних 
алгоритмів при вирішенні будь-яких 
завдань оптимізації та прийняття рішень, 
тому його було обрано для настройки 
синаптичних ваг нейромережі |131. 
Спочатку в запропонованому методі 


відбувається ініціалізація початкової 
популяції хромосомами, що містять 
інформацію | про значення вагових 


коефіцієнтів мережі заданої структури. 
Хромосома складається з генів, які містять 
значення ваг 1 зміщень всіх нейронів 
створеної мережі. У  розроблюваному 
методі використовується дійсне кодування 
для | представлення | значень вагових 
коефіцієнтів у хромосомах. Найкращим 
способом кодування ваг при даному 
поданні топології | є | кодування за 
допомогою матриці, розмір якої 
еквівалентний розміру матриці топології 
ШНМ. Елементи цієї матриці - вагові 
коефіцієнти відповідного зв'язку. Довжина 
хромосоми обчислюється за формулою (2): 


ха0(ТУ)нУ ФО, 1), 0) 


де 0 - кількість нейронів на І-ому шарі; 
Т - кількість ознак у навчаючій вибірці; 
І, - кількість шарів нейромережі | 131. 
Після | формування | проводиться 
оцінювання хромосом поточної популяції, 
для | цього | спочатку | здійснюється 
декодування кожної хромосоми популяції в 
набір вагових коефіцієнтів нейронної 
мережі. Далі виконується обчислення 
значення фітнес-функції, яка оцінює якість 
обраної архітектури за значенням помилки 
навчання нейронної мережі за формулами 


(3) та (4): 


Р -тіл(Х)), (3) 
ко. УЧ9,-в), 4) 
сб і і) о 
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де Сі - потужність множини навчальних 
пар; Ор - отримане за допомогою мережі 
значення вихідного нейрона при і-му 
наборі навчання; К; - необхідне значення 


вихідного - нейрона | при  і-му наборі 
навчання |14). 
Як бачимо, фітнес-функція буде 


оцінювати помилку як різницю між 
отриманим виходом мережі й необхідним, 
яка визначає, наскільки отриманий вихід 
мережі відрізняється від потрібного. Отже, 
чим менше ця помилка, тим вище значення 
придатності, тобто необхідно привести 
матрицю до такого виду, при якому 
помилка навчання нейронної мережі буде 
мінімальною. 

Виходячи | з отриманих значень 
фітнес-функції, відбувається вибір особин 
для генерації нових рішень за допомогою 
рангової селекції. Для цього поточну 
популяцію було відсортовано, відповідно 
до | значень | фітнес-функції, кожній 
хромосомі присвоюється число, що 
визначає її місце в ранзі відсортованої 
популяції, а не абсолютне її значення. 
Даний підхід дозволяє контролювати 
селективний тиск і обмежити кількість 
нащадків однієї хромосоми. Коефіцієнт 
тиску селекції |151. 

Лінійне ранжування усіх хромосом 
виконується за допомогою формули (35): 


Ро5-1 


Р, (Роз) «2-5Р--2.(5Р -1): :0) 


де Ро5 - позиція хромосоми в популяції 
(Ро5-1 має хромосома з найменшим 
значенням фітнес-функції Еіп, а Роз-С має 
хромосома з найбільшим значенням фітнес- 
функції Киах); С - кількість хромосом 
(особин) в популяції; 5Р - коефіцієнт тиску 
селекції, який розраховується за формулою 


(6): 
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де |Епах  -  Хромосома з найбільшим 
значенням фітнес-функції, Бау» - середнє 
значення фітнес-функції всієї популяції 
161. 

Перевага рангового методу полягає в 
можливості його застосування як для 
максимізації, так 1 для мінімізації функції. 
Він також не вимагає масштабування через 
проблеми передчасної збіжності, що є 
актуальною для методу рулетки. 

Оператор схрещування в 
запропонованому генетичному методі 
породжує двох нащадків 0 від двох 
батьківських хромосом, тобто 3 двох 
векторів | дійсних чисел здійснюється 
формування двох нових векторів. У якості 
основи використовується 5ВХ - кросовер, 
який імітує роботу двійкового оператора 
схрещування. Нехай Р'- (р'1,р'2,...,р')) та 
РА-(р'1,р'2,...грих) хромосоми двох предків, 
тоді гени хромосом-нащадків знаходяться 


за формулами (7) та (8): 
1 1 1 2 
су35(1-9)-р) (І нодр)), (7) 
2 1 1 2 
са ко):р) (І-о)рі), | (8) 
де і-1,2, ...,5 2. - довжина хромосоми; ру, 
рі - гени першого і другого батьків, 
відповідно; «о - число, отримане за 
формулою (9): 
(2») у «0,5 
ри 1 
що-то | У», 
21- у) 


де у - випадкове число, розподілене по 
рівномірному закону ує (0, 1); Б - параметр, 
що впливає на ймовірність появи нащадків 
вдалині від батьківських особин Ьє|2,51. 
Під час дослідження будо встановлено, що 
невеликі значення Б дозволяють генерувати 
нащадків, значно віддалених від батьків, 
тоді як збільшення Б спричиняє появу 
нащадків поблизу батьківських пар |17|. 
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У еволюційних алгоритмах оператор 
мутації відповідає за широкий пошук і 
захоплення нових областей пошукового 
простору. Для його налаштування 
користувачеві | додатково | доводиться 
приймати рішення про те, яке значення 
ймовірності мутації слід використовувати. 
При низькій мутації алгоритм буде рідко 


захоплювати нові області пошукового 
простору Й швидше сходитися до 
локальних екстремумів. При високій 


мутації алгоритм частіше обстежуватиме 
різні області пошукового простору, але при 
цьому не буде локалізувати перспективні 
області простору пошуку. Таким чином, 
вибір оператора мутації безпосередньо 
впливає на якість роботи методу 1 
одержувані ним рішення. 

Для вирішення цієї проблеми був 
реалізований прийом зміни імовірності 
мутації. Спочатку на кожній ітерації 
знаходиться рівень розвитку популяції, 
тобто середнє значення фітнес-функції всіх 
особин, після чого отримане значення 
порівнюють із значенням | функції 
придатності кожної особини. Якщо 
значення фітнес-функції індивіда менше 
середнього значення фітнес-функції всіх 
особин, то здійснюється скачок імовірності 
мутації, тобто, чим гірше значення фітнес- 
функції, тим вище імовірність мутації. 
Імовірність мутації хромосоми 
розраховується за формулою (10): 


ові НО) у вир(ду» Р(с), (10) 
г) 


0.5, Е(Р(е))« БК 


Році 


де Епахе) - найбільше значення фітнес- 
функції поточної популяції, К(е) - середнє 
значення фітнес-функції поточної 
популяції (18). 

Якщо ген р; піддається мутації, то 
його | нове | змінене значення рі" 
визначається за формулою (11): 


є 
пеу/ е 
Рі зро, -р,) здо , (1) 


тах 
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де Б - випадкове згенероване число з 
інтервалу |0,1|; О; випадковим чином 
генерується з множини (ріп, ртах), Де Ріпіп, 
і рах - ниЖНЯ і верхня межа можливої зміни 
значення змінної рі; е - номер поточного 
ПОКОЛІННЯ; шах - Максимальне ЧИСЛО 
поколінь; 5 - параметр уточнення, який 
залежить | від | характеру збіжності 
ітераційного процесу після досягнення 
стаціонарного стану, коли кращий з 
індивідуумів | не змінився | протягом 
останніх е поколінь, величина зменшується 
вдвічі, що призводить до розширення 
області пошуку 1, тим самим, до подолання 
пасток локальних екстремумів | 191. 

Така адаптивна мутація дозволяє 
дотримуватися в процесі реалізації ГА 
(еволюції) необхідного балансу між двома 
різномасштабними змінами (мутаціями) 
генів, оскільки на первинних кроках методу 
в основному переважали великомасштабні 
зміни (що забезпечують широку область 
пошуку), в той час як на заключному етапі 
(за рахунок зменшення масштабу мутацій) 
відбувалося уточнення рішення. 

Застосування адаптивної мутації для 
кожного індивідууму є необхідним, але не є 
достатнім для запобігання збіжності 
популяції | в локальному  оптимумі. 
Допустимо, що і-та особина має імовірність 
мутації, яка становить 9790. З одного боку, 
настільки висока мутація при породженні 
від 1-ої особини популяції рівносильна 
породженню хромосоми з повністю 
випадковими параметрами, які не несуть в 
собі корисної інформації про попередню 
еволюцію. З іншого боку, і-та особина 
залишається в популяції через високе 
значення функції пристосованості, 
витісняючи наступні покоління з меншими 
значеннями. Для того, щоб повністю 
уникнути збіжності, введемо селекцію, яка 
буде проріджувати популяцію і відсівати 
особин з високою імовірністю мутації за 
формулою (12): 


їй Виш (Р! ) »бігетоува Р', (12) 
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де Віиш - імовірність мутації Р! хромосоми; 
0 - поріг селекції, що визначає ступінь 
мутацій та відсіює особин, «несумісних з 
життям» |20). 

Метод, що розробляється, зупиняє 
свою роботу при досягненні максимального 
числа епох функціонування, яке задається 
користувачем. 

Представлений генетичний метод 
здатний повністю досліджувати простір 
пошуку, уникати локальних екстремумів на 
етапі генетичного пошуку і ефективно 
використовувати | знайдені «хороші» 
рішення, тобто поступово покращувати 
результати на основі проміжних рішень. 


Результати роботи генетичного 

методу для оптимізації процесу 

десульфуризації димових газів від 
діоксиду сірки 

Для побудови моделі прогнозування 
ефективності процесу знесірчення димових 
газів використано експериментальні дані 
роботи вугільного котла потужністю 1000 
МВТ з 07:00 24 травня 2020 р. по 06:00 31 
грудня 2020 р. Набір складається з 11 
атрибутів та 5330 екземплярів: 

- вміст 80» на вході (мг/м/); 

- температура перед абсорбером (С); 

- температура після абсорбера (С); 

-  об'ємгазу на 1-у сходинку (м//ч); 

- об'єм газу на 2-у сходинку (м//ч); 

- о запиленість після  2-ї сходинки 
(мг/нм?); 

- об'єм газу на рециркуляцію (м//ч); 

- витрата води в абсорбері (м//год); 

-. витрата вапна (т/год); 

-. витрата вторинного реагенту (т/год); 

- вміст 80» на виході (мг/м?). 

Для вирішення поставленої задачі 
було обрано програмне середовище на 
основі мови програмування Руфоп та для 
зручної роботи з масивами даних |і 
формування датасетів використано 
бібліотеку МитРу - пакет Руфоп для 
наукових обчислень. Для будування 
моделей нейронної мережі та роботи з ними 
обрано бібліотеку Кегаз (21) і бібліотеку 
Треапо 221. 
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Важливою вимогою до даних для 
моделювання є їх якість. Якщо дані містять 
так званий "шум", сезонну компоненту, 
викиди, пропуски - це негативно впливає 
на точність прогнозів та якість моделей. 
Також дані, призначені для використання в 
якості навчальних датасетів для ПІНМ, 
повинні бути нормалізовані для зменшення 
похибки та покращення якості тренування. 

Обробка сирих даних перед поданням 
моделі проходить у наступній послідов- 
ності |231: 

-. очищення датасетів з невизначеними 
або пустими ключовими полями; 

- обробка | пропусків в даних 
предикторів; 

-. обробка аномалій даних предикторів; 

- видалення сезонної компоненти з 
часових рядів; 

- приведення даних до типів, 
використовуються у розрахунках; 

-- нормалізація даних. 

При обробці пропусків даних порожні 
значення замінюються на  медіанне 
значення, розраховане за формулою (13): 


що 


р 
ши-3У 
- бала СОНЕНЬ. Др 
М, - Х ме я Ї ( ) 
Ме 
де Хме - НИЖНЄ значення медіанного 


інтервалу; ім - медіанний інтервал; Уме- 
сума спостережень, що була накопичена до 
початку  медіанного інтервалу;  /иг- 
кількість спостережень в  медіанному 
інтервалі. Таким чином, забезпечується 
мінімальна статистична похибка від 
значень ряду. Обробка аномалій в даних 
являє собою очищення набору змінних від 
аномально високих або низьких значень. 
Очищення від сезонної | компоненти 
часового ряду відбувається за допомогою 
методу декомпозиції |23). 

Нормалізоване значення змінної х 
розраховується за формулою (14): 


х - тип (х) 


«х)а 


(14) 


тах (х)- тип (х) 
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Створена мережа проходила навчання 
впродовж 100 епох. При тестуванні було 
використано спосіб поділу вибірки даних 
на навчальну та тестову вибірки у 
відсотковому відношенні 75/2590. 

Для оцінки якості прогнозованих 
моделей використовується середня 
абсолютна помилка (Меап АБзяоїціє Нггог 
(МАЕ)) та середньоквадратична помилка 
(Меап 5дцаге Епгог (М5Б)) (231. 

На рис. І та на рис. 2 представлено 
результати роботи повнозв'язного 
персептрона з одним прихованим шаром. 


тодеї! теап абзоїйіе еггог 


-т бай 
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і 
о 
і 
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Рис. 1. Значення метрики мережі (МАЕ) для 
повнозв'язного персептрона з одним прихованим 
шаром 


З рис. 1 можна побачити, що 
найкраще значення МАН даної моделі 
склало 59,35 та з 10 епохи спостерігається 
практично стале значення протягом 
навчання повнозв'язного персептрона з 
одним прихованим шаром. 
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Рис. 2. Значення метрики мережі (М5Е) для 
повнозв'язного персептрона з одним прихованим 
шаром 
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З рис. 2 можна побачити, що спостері- 
гається поступове зниження значення М5Е, 
а на 20 епосі воно практично стає сталим 
значенням та дорівнює 6785. 

На рис. 3 та на рис. 4 представлено 
результати роботи повнозв'язного 
персептрона з двома прихованими шарами. 
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Рис. 3. Значення метрики мережі (МАЕ) для 
повнозв'язного персептрона з двома прихованими 
шарами 


З рис. З можна побачити, що значення 
МАЄ даної моделі склало 40,98 та, 
починаючи з 20 епохи, спостерігається 
практично стале значення протягом 
навчання повнозв'язного персептрона з 
двома прихованими шарами. 


тодае! І055 


6300) --- їаїп 
того Ві 
6200 1 
6100 4 


6000 1 


055 


5900 1 


5800 1 


5700 1 


т т т т т 
0 20 40 60 80 100 
еросі 


Рис. 4. Значення метрики мережі (М5Е) для 
повнозв'язного персептрона з двома прихованими 
шарами 


З рис. 4 можна побачити, що 
спостерігається | поступове | зниження 
значення помилки, а в районі 10-15 епохи 
тренування досягає локального екстремуму 
і дорівнює 5795. 
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Про локальність мінімуму помилки може 
свідчити подальше поступове зменшення 
помилки мережі. 

На рис. 5 та на рис. 6 представлено 
результати роботи повнозв'язного 
персептрона з двома прихованими шарами 
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Рис. 5. Значення метрики мережі (МАЕ) для 
повнозв'язного персептрона з двома прихованими 
шарами та генетичним методом з адаптивною 
мутацією 


З рис. 5 можна побачити, що значення 
МАЕ для повнозв'язного персептрона з 
двома прихованими шарами та генетичним 
методом з адаптивною мутацією дорівнює 
23,95 та, починаючи з 5 епохи, практично 
стале значення, тому можна вважати, що 
під час тренування досягнуто глобального 
мінімуму помилки, а мережа вважається 
тренованою. 
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Рис. 6. Значення метрики мережі (М5Е) для 
повнозв'язного персептрона з двома прихованими 
шарами та генетичним методом з адаптивною 
мутацією 
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З рис. 6 можна побачити, що спостері- 
гається поступове зниження значення 
помилки, а на 20 епосі воно практично стає 
сталим значення та дорівнює 3253. 

Проведено | порівняльний аналіз 
результатів роботи розроблених моделей з 
такими методами, як: лінійна регресія, 
поліноміальна регресія, логістична 
регресія, метод | найближчого сусіда, 
випадковий ліс, метод мурашиної колонії. 
Критеріями оцінювання були метрики 
МАЕ (Меап АрБвоїше Епог), М5Е (Меап 
Замагей Етгог) та часу виконання (табл. 1). 


Таблиця 1. Порівняльний аналіз результатів 
розробленого методу прогнозування ефективності 
десульфуризації димових газів 


Метод МАЕ | М5Е | Час,с 
Лінійна регресія 79.24 | 9103 363 
Поліноміальна регресія | 77.84 | 9031 401 
Логістична регресія 78.56 | 9041 374 
Метод найближчого 75.76 | 7421 428 
сусіда 
Випадковий ліс 74.15 | 7201 469 
Метод мурашиної | пазо | ягода | 303 
колонії 
Багатошаровий 
перцептрон з одним | 59.85 | 6785 864 
прихованим шаром 
Багатошаровий 
перцептрон з двома | 40.98 | 5795 964 
прихованими шарами 
Багатошаровий 
ЧАРАОНОРОЮ. о ЗАОЮНРЕ ароао арці - о 
прихованими шарами та 
розробленим ГА 


З табл. 1 можна побачити, що най- 


кращі результати за точністю 
прогнозування були отримані У 
багатошарового  перцептрона з двома 
прихованими шарами, його середня 
абсолютна помилка склала 24,95, а 


середньоквадратична помилка дорівнює 
3253. Найгірші результаті були отримані у 
методі з лінійною регресією, її середня 
абсолютна помилка становить 79,24, а 


середньоквадратична помилка дорівнює 
9103. 
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Обговорення результатів розробки 

генетичного методу для вирішення 

завдання оптимізації процесу 

десульфуризації димових газів від 

діоксиду сірки 

Аналізуючи отримані результати 
експериментів (табл. 1), можна побачити, 
що застосування нейронних мереж у 
порівнянні з регресійними моделями дає 
виграш у точності, так наприклад, значення 
середньої | абсолютної | помилки у 
багатошарового  перцептрона з одним 
прихованим шаром або двома прихованими 
шарами відповідно склало 59,85 та 40,98, а 
значення середньої абсолютної помилки 
прогнозування лінійної | регресії | та 
логістичної регресії відповідно дорівнює 
79,24 та 78,56. Отже, нейронні мережі 
можуть досягати більш високої точності 
прогнозування, ніж регресійні моделі, 
оскільки вони набагато краще здатні 
обробляти нелінійну поведінку, але 
нейронні мережі поступаються у часі 
виконанні. Час виконання багатошарового 
перцептрона з одним прихованим шаром 
складає 864с, а поліноміальної регресії - 
дорівнює  401с. Для подолання цієї 
проблеми було розроблено модифікований 
генетичний метод з адаптивною мутацією, 
який | дозволяє підвищити точність 
прогнозування та зменшити час навчання 
нейронної мережі. Під час тестування він 
показав найменше шзначення середньої 
квадратичної помилки прогнозування, вона 
склала 3253, це на 5495 менше, ніж у 
метода випадкового лісу та на 4390 менше, 
ніж у багатошарового перцептрона з двома 
прихованими | шарами. Також час 
виконання багатошарового перцептрона з 
прихованим шаром та розробленим ГА 
дорівнює 564с, що є на 3495 менше, ніж у 
багатошарового  перцептрона з одним 
прихованим шаром та на 4195 менше, ніж у 
багатошарового  перцептрона з двома 
прихованими шарами. 

На підставі проведеної роботи можна 
зробити висновок, що для навчання ШНМ є 
перспективним використання еволюційних 
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алгоритмів, оскільки їх застосування 
дозволяє скоротити час навчання 1 досягти 
більш глибокого мінімуму помилки 
навчання ШНМ. 


Висновки 

У роботі розглядалося вирішення 
проблеми оптимізації процесу 
десульфуризації димових газів від діоксиду 
сірки | з | використанням евристичних 
методів. Було розроблено модифікований 
генетичний метод прогнозування 
ефективності знесірчення димових газів на 
основі | створеної нейронної мережі. 
Принципова відмінність запропонованого 
генетичного методу від існуючих аналогів 
полягає у використанні адаптивної мутації, 
яка використовує в процесі своєї роботи 
рівень розвитку популяції. У ній найменш 
значущі гени з більшою імовірністю будуть 
мутувати в хромосомах, ніж з високою 
придатністю, це підвищує точність та 
збільшує їх роль в пошуку. Тобто, чим 
гірше пристосована особина, тим вона далі 
від оптимуму, отже, Їй надається 
підвищення імовірності мутації, щоб вона 
могла відійти від поточного  неопти- 
мального | положення. Представлений 
генетичний метод здатний повністю 
досліджувати простір пошуку, уникати 
локальні екстремуми на етапі генетичного 
пошуку і ефективно використовувати 
знайдені | «хороші» рішення, тобто 
поступово покращувати результати на 
основі проміжних рішень. Також він 
дозволяє відображати вплив характеристик 
часового ряду параметрів процесу на 
фактичну ефективність десульфуризації та 
підвищувати її точність прогнозу. Було 
проведено порівняльній аналіз роботи 
запропонованого методу та відомих, в ході 
чого було встановлено, що розроблений 
метод дає найменшу помилку 
прогнозування вмісту 502» на виході після 
десульфуризації димових газів (середня 
абсолютна помилка складає 24,95 та 
середня квадратична помилка дорівнює 
3253, що є менше у порівнянні з моделями, 
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побудованими за допомогою інших 
методів) та дає найменший час виконання 
(час навчання моделі склав 564 с). 
Практичне використання розробленого 
методу дозволить підвищити ефективність 
процесу десульфуризації димових газів та 
зниження викидів 502 в атмосферу. 
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